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STRUCTURAL  OF  FLAME  H;>  ~  Oj  UNDER  PRESSURE 


FINAL  REPORT 


The  contract  between  the  ERO  London  and  Padua  University  -  Institute 
of  Chemical  Plants  refers  to  "Studies  of  the  structure  and  composition  of 
a  flat  hydrogen-oxygen  flame  under  pressure  by  means  of  thermocouples  and 
spectroscopic  methods.  The  experimental  results  will  be  raathema tically 
correlated" . 


ABSTRACT 


The  report  illustrates  the  work  carried  out  over  5  years,  summarizing 
and  correlating  that  which  has  already  been  set  out  in  eight  technical 
reports  on  the  state  of  progress  of  the  research  with  the  addition  of  how 
much  has  been  done  since. 

Particularly  described  are:  the  reactor  and  subsidiary  equipment , the 
equipment  and  measuring  methods  of  temperature  profiles,  concentration  of 
radical  species,  combustion  velocity  and  finally  a  new  hybrid-stochastic 
method  of  flame  simulation. 

Particular  emphasis  is  given  to  how  much  has  newly  been  done  from  mi 
experimental  and  theoretical  point  of  view  and  to  the  results  obtained  in 
a  field  "flame  under  pressure",  heretofore  studied  little, and  in  which  the 
results  given  in  literature  are  often  contradictory. 
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a )  Reactor  and  subsidiary  equipment. 

The  lay-out  of  the  reactor  built  by  us  and  used  in  the  study  of  flames 
under  pressure  is  given  in  Fig.l. 

The  characteristic  parts  are: 

-  the  burner  with  corrugated  matrix;  which  ensures  a  perfect  mixing  of  the 
entering  mixture  and  laminar  plug  flow  in  the  reactant  mixture. 

-  the  combustion  chamber,  =  5.6  cm,  a  double  transparent  quartz  tube 

that  allows  working  pressure  uf  up  to  6-7  10^  Pa  and  the  carrying  out  of 
optical  and  spectroscopic  measurements  and  observations  of  flame  geometry 
and  on  the  position  of  probes  and  thermocouples. 

-  the  flameholder  is  a  holed  double  disk  in  steel  and  pyrophy lli te . 

The  combustion  chamber  is  equipped  with  automatic  scanning,  Fig. 2  wich 
allows  on  one  hand  a  high  spatial  resolution  in  the  measurements  by  means  of 
a  thermocouple  introduced  in  the  flame  as  well  as  in  the  optical  and  spectre 
scopic  ones, and  on  the  other  velocity  in  the  experimental  surveys  thus 
ensuring  constancy  in  the  behaviour  of  the  flame. 

Fig.  3  shows  the  arrangement  of  the  burner  and  instruments  both  for 
combustion  speed,  flame  temperature  and  concentration  of  chemical  species. 

Fig.  4  and  5,  in  particular,  show  in  detail  the  gas  flow  system  and 
the  optical  instrumentation. 

The  various  pieces  of  equipment  are  described  in  the  paragraphs  that 
concern  the  control  of  the  progress  of  the  reaction  or  measurement:  particulars. 

b )  F ] o w  Velocity  and  P  r  e  s s  u  r  o  Me a  s  n r o m o n  t  s . 

Next  to  tlie  reactor  a  critical  orifice  is  set  up  which  produces  a 
pressure  drop,  so  that  the  pressure  in  the  combustion  chamber  is  maintained 
at  a  constant  value.  Critical  orifices  are  also  placed  in  the  feeding,  lines, 
these  allow  the  flow  velocity  of  the  single  components  and  the  composition 
of  tlie  mixture*  to  the  burner  to  lie  measured.  The  pressure  before  and  aftei 
is  measured  by  transducers  of  the  strain-gauge  typo  (Philips  Type  PR  9F-0>) 
set  up  in  Whe.it  sl«ne  bridges  on  a  thim  sheet  of  steel. 

The  sig.ua  Is  of  the  tin  re  feeding,  lines  and  that  at  tin*  m::v*r  millet 
are  measured  on  a  iiuriiift  ir  "Tape"  simultaneously  with  the  t  rmpe  ra  i  mv  and 
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pressure  values  in  the  combustion  chamber.  The  signals  are  also  visualised 
on  milliameters  set  up  on  the  operative  bench. 

The  transducers  can  visualise  variations  of  few  ram  of  Hg  in  normal 
working  conditions  up  to  10.10^  Pa. 

c )  Combustion  velocity. 

It  is  known  how  this  parameter  is  one  of  the  most  important  in  flame 
study.  Data  in  the  literature  are  scarce,  in  particular  for  combustion  under 
pressure.  Experiments,  if  conducted  with  accuracy,  are  the  most  trustworthy 
to  obtain  data  on  the  combustion  velocity  of  the  mixture  under  examination. 

The  measurements  are  conducted  visualising  the  flow  by  means  of 
suspended  particles  of  MgO.  The  stroboscopically  illuminated  particles  arc 
photographed  every  two  seconds  with  0.5  seconds  exposure  time.  The  camera  is 
an  Asahi-Pentax  with  a  135  mm  P  2.5  telephoto  lens.  A  double  enlargement  is 
obtained  and  a  focal  range  from  1  to  2  rum.  The  range  is  focused  on  the  luminous 
particles  that  are  found  on  the  axial  diameter  of  the  flame. 

The  velocity  is  determined  by  diagramming  the  height  of  the  particles 
against  time.  The  tangent  to  the  curve  at  a  point  gives  the  combustion 
velocity  at  that  point. 

d )  Temperature  Profiles. 

The  temperature  profile  of  a  flame  is  that  which  analyses  more  completely 
its  evolution,  visualizing  the  preheating,  reaction  and  post-reaction  zones. 

The  temperature  values  along  the  flame  are  obtained  through  measurements 
by  means  of  thermometric  probes  immersed  in  the  flame  and  of  special  dimensions 
and  shape,  in  order  to  disturb  the  gas  flow  as  little  as  possible. 

The  probe  used  is  a  25  microns  diameter  chrome 1 -a lumel  thermocouple 
covered  with  a  fused  mixture  of  borax  and  alumina  to  eliminate  the  catalytic 
effects  as  much  as  possible.  The  very  small  diameter  and  t.he  automatic 
scanning  allow  for  good  spatial  resolution. 

The  instant  position  of  t lie*  probe4  is  recorded  by  the  aid  of  a  moving, 
transducer  which  is  capable  of  resolving  the  position  with  .H  30  microns 
accuracy.  The  effective  measurement  time  is  a  few  seconds, so  as  to  ensuvi 
the  constancy  of  the  operative  variables.  Owing  to  the  brief  nva sui. emeu t  t  *, 
it  is  not  possible  to  use  a  normal  pen  iveordor,  whose  constant  is  high,  r 
than  that  of  t  hr  mtunircii'oji  I  ,  ihtoeforr  a  magnetic  tape  record,  r  (Ana  J  op,  7 
Philip'.)  is  used.  Tin*  millivolts  are  t  i  a  sn  1  pi  i:ud  through  a  polynomial  oi 
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the  7th  order  which  fits  the  data  of  the  readings  and  the  temperature 
obtained  is  then  corrected  for  radiation  losses  of  the  thermocouple  junction. 

e )  Concentration  measurements  of  the  OH  and  H  species. 

The  concentration  values  of  the  radical  species  are  an  indispensable 
complement  to  the  temperature  measurement  and  a  useful  indication  on  flame 
behaviour,  in  particular  in  the  reaction  zone  where  the  pressure  effects  can 
have  a  great  effect  on  the  reaction  kinetics. 

The  relatively  low  flame  temperatures  c^°  not  a^ov,T  t^ie  aPpi^-“ 

cation  of  the  emission  spectroscopy  except  in  the  case  of  the  fundamental 
lines  of  alkaline  metals  injected  in  the  flame.  This  technique  permits  the 
determination  of  the  H  concentrations.  The  emission  intensity  of  the  sodium 
D  lines  is  due  only  to  the  chemiluminescence  according  to  the  reactions: 

H  +  H  +  Na  - H9  +  Na 

H  +  OH  +  Na - P*~  H20  +  Na 

since  the  thermal  emission  is  irrelevant,  because  of  the  low  temperatures 
of  the  flame. 

On  the  other  hand,  absorption  spectroscopy  is  used  for  OH  concentration 
measurements.  For  such  measurements  a  quartz  xenon  lamp  is  used  as  light 
source  and  a  Philips  XP  1003  photo-multiplier  as  transducer. 

To  raise  the  spatial  resolution,  cylindrical  lenses  focusing  on  the 
centre  of  the  flat  flame  arc  used.  The  image  on  the  entrance  of  the  spectrometer 
a  (Jarrel-Ash)  is  enlarged  twice.  The  Oil  spectrum  band  (0,0)  from  3060  to 
3200  A,  is  taken  in  the  2nd  order  to  raise  the  resolving  power  of  the 
spectrometer. 
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2.  -  MATHEMATICAL  TREATMENT. 


a )  Heat  Flow. 


The  volumetric  chemical  heat  release  rate  results  at  any  point 


d  dt 

-r-  (m  c  T  -  X  — ) 
dz  p  dz 


(1) 


m  =  p  v  -  mass  flow  rate  of  the  mixture  (p  ,  v  density  and  burning 
o  o  o  o  ° 

velocity  at  normal  conditions) 

0^,  X  =  specific  heat  and  local  thermal  conductivity. 

Equations  (1)  describes  the  heat  release  rate  as  the  gradient  of 
energy  flow  due  to  the  mass  transfer  (convection)  and  to  the  heat  conduction. 
It  does  not  make  into  account  the  heat  losses  which  are  due,  in  particular, 
to  the  radiation.  The  equation  has  been  numerically  calculated  by  Stirling's 
method. 


The  integrated  heat  rate  througouht  the  flame  interval 

Z  =  oo 


%  " 


q  dz 


(2) 


z  =  0 


can  be  compared  (apart  from  the  small  heat  losses)  with  the  heat  calculated 
from  the  reaction  enthalpy  A 

%  ■  <Vvtot>  ”o  vo  <3> 

V  »  V  =  volumetric  flow  rate  of  0o  and  total, 
o  tot  2 

A  comparison  between  the  two  heats  allows  good  agreement  with  the 

temperature  measurements. 


b)  Stochastic  hybrid  model  of  flame. 


In  the  mathematical  simulation  of  the.  burning  process  of  a 
monodi.nonsional  flame  stochastic  solution  lias  boon  adopted  for  the  chcinicaJ 
process,  which  has,  a  more  consistent*  physical  basis  than  the  commonly 
used  do  termini  slic.  It  adapts  better  to  the  use  of  the  computer  (avoimlJir- 
11  round -of f "  errors)  and  allows  the  use  of  calculating  means  to  render  the 
adopted  model  more  exact  and  concise. 

The  storaslir  solution  presents  a  gnod  elasticity  so  much  «:o  Mi.  1  it 
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has  been  possible  to  hybridise  it  with  deterministic  solutions  for  energy 
transfer  and  molecular  diffusion. 


Volume  V,  spatially  homogeneous  and  in  thermal  equilibrium,  contains 
a  mixture  of  various  chemical  species  which  react  among  themselves  in 
different  ways.  They  are: 


N.  =  number  of  particles  of  the  various  chemical  species  S^: 


i  1,2, . ,n. 

Rf  =  chemical  reactions,  r  =  1,2, . ,ra. 


The  probability  of  the  reaction  taking  place  is: 

Vr>  =  V  2  Qr 

i 

and  the  probable  time  is: 


(4) 


m  m 

P  U)  =  £  Q_  exp  -  (2.  Q  t)  (5) 

1  r  l  r 


where  Q  =  c  f  ,  with  f  combinatory  function  of  N.  and  c  ,  stochastic 
r  r  r  r  -  J  i  r 

•  p—  1 

reaction  constant  =  q!  ,  q  =  1,2,3  for  differents,  two  or  three 

equal  reacting  molecules  and  g  =  reaction  order. 

In  order  to  select  the  event  and  the  probable  time  it  is  necessary 
to  generate  two  casual  numbers  and  included  between  0  and  1  and  to 
choose  r  so  that: 

r-1  m  r_ 

Z  i  <6> 

i  i  i 


and 


t  =  0/  Z  Q  )  In  (1/L  ) 
1  i 


(7) 


Once  r  is  known  the  values  of  1L  are  updated  at  time  t  according  to  the 
s t o i c h i ome try  of  rea c t i o n  R  . 

Once  the  reaction  enthalpies  and  the  heat  capacities  of  the  various 
components  are  known,  we  have  for  every  calculation  cycle  s,  and  for  the 
adiabatic  flame: 

T  =  T  -  (AH)  />_  (c  N./N)  k  -  0,1,2,'.... 

b  b  J.  •  •  pi 


and  l  lu*  average  flow  velocity  will  be: 

(NT) 


V  -  V 


s-J  (Nlj 


O 
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The  conservation  equation  for  the  species  i  in  the  case  of  "premixed, 
stationary,  laminar,  ideally  monodimensional  flame"  is: 


d 

dz 


<v+V 


■] 


dn. 
_ 1 

dr 


(8) 


with  z  =  axial  coordinate,  n  =  mole  density  number,  v  and 
=  mean  velocity  and  diffusion. 


Introducing  the  fractional  mass  flux  we  have: 


P  v 
o  o 


M. 

l 


d  G. 

l 


=  d  n. 

i 


(9) 


with  p  v  =*  mass  flux  velocity 
o  o  J 

=  weight  fraction  of  species  i  in  the  flux. 

M  =  molecular  weight. 

Since  the  variation  of  takes  place  in  a  stochastic  way,  per 
unitary  increase,  the  displacement  of  one  molecule  will  happen  when: 


p  V  G.  (  >1  negative  diffusion 

-2—^  -i  |  (10) 

fb  v  v  <  -1  positive  diffusion 

For  the  type  of  flame  adopted  the  energy  conservation  equation, 
integrated  with  the  boundary  conditions  z  =  0,  T  =  Tq,  3T/3z  **  0  offers 


AT 

Az 


v  p  c  (T  -  T  ) 
o  o  p  o 


(11) 


with  0^  and  A  -  mean  specific  heat  and  average  conductivity  of  the  mixture. 

Equation  (11)  is  called  upon  only  for  unitary  events  linked  with  the 
initial  reactions  between  molecules  (reactions  9  and  10  of  Table  2) .  When 
these  reactions  which  arc  the  most  probable  in  the  high  part  of  the  flame 
arc  defined  by  the  programme,  their  stochastic  time  determines  the  entity 
of  heat  transfer. 

In  this  way  both  for  molecular  diffusion  as  well  as  heat  transfer  a 
system  of  hybridization  is  chosen  adopting  the?  deterministic  equations 
only  for  unitary  events  linked  to  the  stochastic  approach. 

Thus  it  is  possible  to  calculate: 

-  the  number  o<"  events; 

-  the  temperature  and  concentration  profiles; 

-  the  slochaslic  time  velocity  and  tin*  distance. 


3. 


RESULTS. 


Among  the  numerous  runs  described  in  the  8  reports,  two  have  been 
chosen  which  characterize  better  the  various  theoretical  and  experimental 
aspects  and  on  the  other  hand  are  sufficient  in  describing  the  behaviour 
of  a  flat  flame  under  pressure. 

The  characteristics  are  given  in  Table  1. 


TABLE  1 


Chamber 

pressure 

Ingress 

Composition 

H?  °2  N2 

Flow 

Velocity 

Qi_ . 

eq(3) 

Qs. 

eq  (2) 

% 

cm  sec 

1 

16.15  5.40  78.45 

9.05 

2.80 

2.55 

3 

12.40  5.15  82.55 

12.85 

11.30 

10.51 

. Characteristics  of  flame  at  atmospheric  pressure  and  under  pressure. 

As  regards  the  heat  release  in  general,  a  discrepancy  of  almost  10% 
is  found  between  the  values  calculated  by  equation  (3)  and  those  deduced 
by  the  temperature  profile  equation  (2).  The  agreement  is  good,  taking 
into  account  the  the  heat  due  to  radiation  is  not  calculated. 

The  reaction  probability,  the  temperature  profiles  and  the  composition 
of  the  two  flames  have  been  calculated  with  the  help  of  a  Laben  70 
computer . 

It  is  assumed  N  =  50.000. 
tot 

Of  the  38  reactions  that  can  happen  in  flames  ^2~^2~^2  on^  ^  are 

necessary  to  obtain  satisfactory  results  from  the  simulation  without  making 

the  calculation  too  involved .The  10  reactions  are  given  in  Table  2  together 

with  the  kinetic  constants.  These  last  are  suggested  by  BAULCH-DUYSDALE  - 
2 

HOME -LLOYD  (  ) . 

Reaction  enthalpies,  specific  heats  and  thermal  conductivities  are 

3 

those  given  in  the  table  by  FENNER  (  )  except  that  for  the  values 

4 

given  in  ref.  (  )  have  been  used. 


Aj 
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TABLE  2 


Reaction 

,  .  ,  3  -1  -1. 

k  cin. (cm  mol  sec  ) 

AT,K 

AH  298(kJ/mole) 

1) 

02+H 

0+OH 

2.2*10U  exp(-8450/T) 

700*2200 

70.672 

2) 

h2+o 

+  H+OH 

1.8«1010  T  exp(-4480/T) 

4 00 t 2000 

8.255 

3) 

h2+oh 

+  E20+H 

2.2*1013  exp(-2590/T) 

300*2500 

-63.304 

4) 

h+ho2 

•*  Oli+OH 

2.5*10^  exp(950/T) 

290*800 

-159.979 

5) 

h2o+h 

+  H2+0H 

9.3*1013  exp(-10250/T) 

300*2500 

63.304 

6) 

1I+H+N, 

r  v»2 

2.5*1018/T 

300*53 00 

-435.973 

7) 

02+H+N2-*H02+N2 

1.5*1015  exp(500/T) 

300*2000 

-197.096 

8) 

h+ho2 

+  H20+0 

5*1012  T°'5  exp(-2000/T) 

300*350 

-231.538 

9) 

V°2 

-*■  OH+OH 

8«1014  exp(-22500/T) 

298*2000 

78. 927 

10) 

V°2 

+  H+H02 

5.5'1013  exp(-29100/T) 

298*2000 

238.906 

The  diffusion  velocity,  owing  to  the  large  excess  of  in  the 
flame,  is  approximated  by: 

ND.  M0  d  N. 

(v  ).  - - _ -  (1 

^  d'i  N.  dz  N  u 

J  *  - 


The  binary  diffusion  coefficient  is  evaluated  from  the  well-known 


equation. 


D.  „  =  2.628 

i.N2 


T3/2  C(Mi"^2)/(2  Mi 


0  *  * 

of.  p  a  (t  ) 

ij 


—  10  ^  cm*"  sec  ^  (13) 


The  values  in  the  constants  are  given  in  ref.  (  )  while  the 
corrections  caused  by  the  presence  of  polar  species  arc  valued  according 

6  jj- 

to  (  ) .  In  the  simulation  programme  the  reduced  collisional  integral  ft 

*  -1  /  2 
(T  )  is  fitted  with  a  10th  degree  polynomial  expression  in  (T)  from 

7  .  .  . 

the  data  in  (  ) ; the  polynomial  coefficients  are  given  m  Table  3 

TABLE  3 

7^0^401/, 9  b  ^  1  .036917  c  »  -0.036973 

d  =  0.818693  e  «  -1.416095  f  -  2.153107 

g  =  6,386901  h  =  -2.910819  i  -  -7.933206 

1  -  -3.845663  m  ^  9. 390/0 J 

The  programme .  are  given  in  the  appendix. 


-  16  - 


a)  Contribution  of  every  single  model  reaction. 

The  calculation  programme  takes  into  account  the  stochastic  path  in 
the  form  of  number  of  event  of  the  various  elementary  reactions  carried 
out  by  the  same  programme  (Fig. 6).  Such  data  have  allowed  the  programme 
optimization  to  reduce  notably  the  troublesome  kinetic  model  discaring 
the  radical  reactions  which  in  our  case  were  practically  nil. 

From  the  diagrams  of  Fig. 6  it  is  possible  to  observe: 

-  the  reactions  between  radicals  tend  to  diminish  at  the  approach  of 
equilibrium.  This  can  be  represented  by  reactions  (3)  and  (5) . 

3 

H2  +  OH  H20  +  H 

5 

that  reach  in  proximity  of  equilibrium,  intervention  values  near  to  100%. 

-  reaction  (9)  between  H2  and  02  tends  towards  zero  corresponding  with 
the  ignitions  point  of  the  flame. 

-  reaction  (4)  shows  a  maximum  at  intermediate  temperatures,  a  zone  in  which 

the  species  H  02  is  stable,  and  which  later  decreases  to  zero  at 
equilibrium.  v 

b )  Temperature  Profiles. 

The  temperature  profiles  are  given  in  Fig. 7.  A  dotted  line  shows 
the  experimental  ones,  a  broken  line  indicates  the  calculated  ones.  The 
agreement  can  be  regarded  as  satisfactory  taking  into  account  that: 

-  the  experimental  values  suffer  from  approximate  measurements  and  an 
imperfect  flat  shape  of  the  flame. 

-  in  the  model  the  flame  is  assumed  adiabatic  and  certain  secondary  reactions 
are  disgrcgarded. 

From  the  two  profiles  of  1  and  3  atm  it  is  possible  to  note,  in 
agreement  with  the  increase  of  the  global  reaction  velocity  with  pressure: 

In  (r^/r2)  =  £(q/2)  -  lj  In  (p^/p2)  where  q  is  the  global  reaction  order  - 
=  2.081  (^),  that  equilibrium  is  reached  more  quickly  as  the  pressure  is 
increased  and  the  reaction  zone  becomes  thinner. 

It  is  also  noted  (runs  described  in  previous  reports),  that  the 
effect  of  diffusion  is  more  accentuated  at  low  pressure  than  at  high 
pressure,  given  that::  tin*  global  reaction  velocity  is  proportion'll  to 
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p(q/2)_i  .  t^e  diffusion  is  inversely  proportional  to  pressure,  as  it 
increases,  the  reactions  gain  in  importance. 

c)  Concentration  Profiles. 

The  facts  previously  described  find  verification  also  in  the 
concentration  profiles,  Fig.  8-9.  The  increase  of  the  global  velocity 
leads  to  a  quicker  reaching  of  equilibrium  conditions. 

The  two  most  important  chain  carriers,  OH  and  H,  show,  as  intermediate 
species,  concentration  maxima  whose  position  depends  on  pressure.  The 
relative  curve  at  OH  presents  two  of  them,  the  first  relative  to  reaction 
(4)  masimum  and  the  second  together  with  that  of  H,  due  to  reactions  (1) 
and  (2)  maximum. 

The  presence  of  a  minimum  for  OH  had  been  noted  in  the  experimental 
determination  of  such  species,  report  6. 


CONCLUSIONS  . 


The  theoretical-experimental  study  of  flames  unc*er  Pressurc 

has  allowed  us  to  identify  the  influence  of  the  pressure  parameter  on 
fluid  dynamics  and  on  the  progress  of  the  system  reactions  in  flow. 

The  small  reaction  zone  of  flame  under  pressure  and  the  low  temperature 
of  flame  due  to  the  excess  of  are  favourable  conditions  in  the  assumption 
of  single  dimensional  flow  and  therefore  ideal  for  the  experimental  chech  of 
the  simulation  model  adopted. 

This, a  hybrid-stochastic  type,  offers  better  opportunities  in  respect 
of  the  deterministic  normally  used.  In  fact  besides  avoinding  "round  off 
errors",  giving  stable  and  accurate  solutions,  it.  is  capable  of  indicating 
the  most  important  physical  phenomena  and  the  most  probable  reaction 
mechanism  for  the  various  temperature  and  pressure  conditions  of  the 
flame  considered. 
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